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^ , Abstract 

. We consider the ensemble of NxN hermitian random matrices whose 
entries are equal to zero outside of the band of width b along the principal 

. diagonal. Inside this band the entries {Hij, i < j} are given by independent 

' identically distributed gaussian random variables with zero mean value and 

. variance /b. 

' We study asymptotic behavior of the spectral norm in the limit 

. N ~^ oo when the band width b is much smaller than the matrix size N 

Qh' but also tends to infinity. Our main result is that if 6/(logiV)'^ oo, then 

(— I . limsup^ J is bounded by 2v with probability 1. To prove this, we de- 

"j^ , rive a system of recurrent relations for the moments M^^'''^ and analize these 

relations in the limit when k goes to infinity as iV ^ oo. 

^ ; 1 Random matrices and eigenvalue distribution 

H , 

I Random matrices of infinitely increasing dimensions N ^ oo were considered first 

in theoretical physics. Physics serves as the source of new questions and results here 
and the development of the random matrix theory reveals many relations between 
random matrices and other branches of mathematics (see the special issue on physics 
and mathematics of random matrix theory of the Journal of Physics A : Mathemat- 
ical and General, vol. 36 (2003)). In these studies, the special attention is paid to 
the extreme eigenvalues oi N x N hermitian random matrices with unitary invariant 
probability distribution. Recent progress in this field relates random matrix theory 
with nonlinear differential equations and integrable systems (the pioneering mathe- 
matical results are obtained in [23); from another hand the probability law of the 



maximal eigenvalue of random hermitian matrices ImaL is asymptotically equivalent 
to that of the longest increasing subsequence in the random permutation of 1, . . . , 
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(see e.g. El)- These results motivate extended studies of extreme values of 
random matrix ensembles. 

In present paper we study asymptotic behavior of /max of the ensemble of band 
random matrices. This ensemble plays an important role in the quantum chaos 
theory (see e.g. EH). The probability distribution of this ensemble is not invariant 
with respect to the unitary transformations. 



1.1 Wigner ensemble and the semicircle law 

The spectral theory of random matrices was started half a century ago by E. Wigner 
(see e.g. |24|'). He studied the eigenvalue distribution of the ensemble oi N x N real 
symmetric random matrices W^^^ of the form 



where {wij,i < j} are independent random variables. The main proposition proved 
in |24j is that if Wij are centered random variables of the variance v"^ , odd moments 
zero and all even moments finite, then the normalized eigenvalue counting function 
of W^^^ defined by the formula 

^7v(0 = < Oiv^S 

where ^ < . . . < are eigenvalues of W^^^ , weakly converges in average to the 
limiting function a{l; 2v) 

lim crAr(/) = a{l;2v) (1.2) 
with the density of the semicircle form; 



aV, 2v) = Pit; 2v)='{ f"' f ■ ^ (1.3) 
Zttv^ I U, otherwise. 

This statement is known as the semicircle (or Wigner) law. In particular, Wigner 
has shown that the moments of the limiting distribution are determined by 
recurrent relations 

k-l 

m2k = v'^Yl "^2fc-2-2i m2j, mo = 1. (1.4) 
jo 

Obviously, m2k+i = 0. 
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1.2 Gaussian Unitary (Invariant) Ensemble 

If one considers the ensemble of N x N hermitian matrices 

= -^h,,, (1.5) 

where hij = Oij + i(3ij,i < j are jointly independent complex gaussian random 
variables with zero mean values and variance v'^, the density of the probability 
distribution of H'-^^ = H can be written in the form 

Z-^e^^i-^TiH^], (1.6) 

where is the normalization constant. Note that in this case random variables a 
and (3 are jointly independent (for more details, see section 5 of the present article). 

The probability distribution (1.6) is invariant with respect to the unitary trans- 
formations. Therefore the random matrix ensemble (1.5)-(1.6) is referred to as the 
Gaussian Unitary Ensemble (GUE) of random matrices. It plays the central role in 
the spectral theory of random matrices (see e.g. j2Uj ) . 

It is easy to see that j//'^^-*} verifies the conditions imposed on {IV^^^-'} and 
therefore the semicircle law (1.2)-(1.3) is valid for GUE. Moreover, convergence 
(1.2) holds with probability 1 in this case. The first study of the maximal eigenvalue 
lraa.yi{H^^^) = \\H'^^'^\\ was Carried out by S. Geman m for the ensemble of real 
symmetric matrices having the probability distribution of the form (1.6). This en- 
semble is known as the Gaussian Orthogonal Ensemble (GOE) of random matrices. 
Geman has proved that /max 

converges with probability 1 as — > oo to the 
border 2v of the support of da{l; 2v). 



2 Band random matrices 

Given a parameter 6, let us consider hermitian random matrices H^^'^^ of the form 



^iv''^ = /^xW^S, x,y = l,...,N, (2.1) 



^xy — "'xy M ^xy 

where {hxy, x < y} are the same as in (1.5) with the law (1.6) and 

'X - y 



^xy t \ 



where 



, , , 1, if |t| <l/2, 
^ '0, otherwise ^ ' 
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The ensemble of band random matrices (2.1)-(2.2) was considered in relation 
to the quantum chaos theory and solid state physics (see [SI 1121 llUj). The most 
intriguing question was that the ratio b'^ /N separates two major asymptotic regimes 
that characterize behavior of certain spectral characteristics of random matrices. 

In paper ^7j, the real symmetric analog of (2.1) is considered. It is proved 
that in the limit b — > 00,6 = o{N) the Wigner law is valid, i.e. the normalized 
eigenvalue counting function aisifi{l) of H^^'^^ weakly converges in probability to the 
semi-circle distribution (see also papers j7|El|^ for this and related result). The 
same proposition can be easily proved for the ensemble of hermitian band random 
matrices. If 6 = 0{N''') with some 7 > 0, then the convergence of <T7v,6 to the 
semicircle distribution holds with probability 1. 

Much less is known about the limiting behavior of the spectral norm 

Imax of band random matrices. Up to our knowledge, this question was addressed 
in paper |5] only. If one looks at the computations presented in [S], one can easily 
see that — > cxd as 6 = o{logN),b 00. In present paper we show that 
if b goes to infinity faster than (logA^)^, then the spectral norm of band random 
matrices remains bounded by 2v. 

Our main result is given by the following proposition. 

Theorem 2.1. 

If N and b go to infinity in the way such that b/(logN)^ — > 00, then the spectral 
norm remains hounded with prohahihty 1; 

limsup||F(^'^)|| < 2w. (2.4) 

Af,b— >oo 

Ifb = OiN'f) with 7 > 0, then (2.4) turns to equality 

To study the spectral norm of random matrices H, we employ the general ap- 
proach used first by S. Geman by suggestion of U. Grenander Later it was 
employed by many authors in applications to various random matrix ensembles with 
jointly independent entries OIHIEI and also for random matrices whose elements 
are statistically dependent random variables 0]. The key observation is that if one 
considers the moments 

M(f)=E{lTr[i/Wp} 

in the limit N 00, k = 0{\ogN), then the leading contribution to Mg^^ will be 
given by the maximal eigenvalue of /max(-ff)- 

The next important step was made in [3]. It was observed that to study the 
limiting behavior of the moments M2'^^ of the Wigner ensemble (1.1), it is not 
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necessary to compute them explicitly. It is sufficient to show that they verify a 
system of recurrent relations that converges as N,k ^ oo to the system of equalities 
(1.4) that determines the moments of the semicircle law. 

In paper we have developed a method of derivation and study the recurrent 
relations for the moments in the case when the matrix elements are gaussian 
correlated random variables. In present paper we develop a new version of this 
method that we adapt to the band random matrix ensemble. Also we present more 
precise analysis of these relations aiming the best possible estimates of the moments 
M^^f') of H(^'b\ 

We start with the basic example of the GUE (see Section 3). We pay much 
attention to this ensemble because the relations and estimates of the moments Mg^'''^ 
are very similar to those of GUE. Thus, the estimates we get in the GUE case are 
true for the moments of the band random matrices. We did not manage to get the 
optimal estimates for Mg^'* obtained by other methods. The benefit of the approach 
developed is that it can be applied for random matrix ensembles different from GUE. 

In Section 4 we modify our approach and study the band random matrices. These 
computations lead to the proof of Theorem 2.1. In Section 5 we prove auxiliary 
statements. 



3 Moments of GUE matrices 



3.1 Recurrent relations 



In this section we derive two main recurrent relations for the moments of H^^'^^ and 
their variances. Let us consider the normalized trace La = jfTrH"' and compute 
the mathematical expectation with respect to the Gaussian measure (1.4) of the 
following expression, where G denotes some regular function of H: 

1 ^ 

x,s=l 

Regarding the last average, we can use the integration by parts formula 

^f— |. (3.1) 

The use of the partial derivative symbol is explained in section 5. 
Then we can write that 

2 N 

x,s=l j=Q 
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a-2 2 N 

v'J2^{La-2-jLjG} + ^ ^ B{H:-'dG/dH,,}. 

j=0 x,s=l 

If G = 1 and a = 2k, then we get the following relation 

2k-2 

j=0 

Denoting = L — E L, we can write that 

E {VL„J = EL„,EL„, + L>g^^ = E{L°^L°J. 

Then using the fact that ^2^^ = EL2fc+i = (see Section 5), we obtain our first 
equality 

= E <i-2.Mif ) E I^g., (3.2) 

j=0 ai+a2=2k-2 

that form the system of recurrent relations. Let us note that for two random variables 
we always have E (LfL^) = E (LiL^). 

Following our general scheme, we introduce variables 

Dit.,a, = E {K^K^ • • • = E {LaAK, ■ ■ ■ Kj} 

and apply (3.1) to this expression with G = [L°^ ■ ■ ■ L°J°. Then we obtain relation 

j=0 

^,2 9 ^ ^ o "1 

— ^E jL^^ • • • Ll._^ Qi -Lai+ai-2 ■ ■ ■ K,\ ■ 

1=2 

The last term arises because of equality 

Now we use two times the identity 

E {L1L2Q} = E {LiL^Q} + E {LIL2Q} + E {LlL^Q} - E {L?L^}E {Q} 

and obtain our second recurrent relation 

ai— 2 ai— 2 

-^ai,...,Og - ^ 2^ -""^j ^ai-2-j,a2,...,ag ^ ^ ^^^ai-2-j^j,a2,...,ag^ 
j=0 j=0 



E ^ = E E = («-!) (w^'^-O, 
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ai— 2 ai— 2 

2 V- _ 2 r)(2) 

^j,ai-2-j,a2,...,a, -^j,ai -2-jr-^a2,...,a, ^ 

j=0 i=0 

/Y2 .iL^ « ai+aj-2-^a2, ai+i,.--,a<7 ^ l\f2 * a2v,ai_i,ai+ai-2,a,+i,...,ag- V'-'-"-'; 

i=2 i=2 

When deriving (3.3), we assumed q > 3. For q = 2 we have equahty 

ai—2 ai— 3 

^Sa2 = E Mr^a,-2-i,a2 + E ^£\-2-,,a2 + ^v'Ma,+a,-2. (3.4) 

i=o j=i 

In the case of g = 3 we adopt relation (3.3) assuming that = 1 and D^^) = 0. 
Also one has to remember that Mg^^^ = 0. 

3.2 Recurrent estimates 

Let us prove the following proposition. 
Lemma 3.1 If N > 4 and k < N^/^, then 

k 



and 



m2k (3.5) 



|I^il..J<^i^(l + ^j -2. (3.6) 
where 2k = ai + . . . + aq. 

We prove this statement following the procedure of recurrent estimates proposed 
in [Ij. Let us consider the plane of positive integers, i.e. the set of points B = 
{{l,p) ■ l,p € N}. Actually, S is a quarter-plane. 

We say that the random variable D^a},...,aq with given Aq = (ai, . . . , aq) belongs 
to the point {2k, q) € -B, if the sum of a's is equal to 2k. We call such variables D 
the elements of {2k, q). Obviously, 2 < q < 1,1 = 2k and we restrict ourselves with 
the corresponding subdomain B of B. 

It follows from relations (3.3) that the elements D of {2k, q) are expressed in 
terms of the elements D belonging to one of the four other points: {2k — 2,q), 
{2k — 2,q + 1), {2k — 2,q — 1) and {2k — 2,q — 2). If any of these points is situated 
outside of the domain B, then the corresponding term in (3.3) is equal to zero and 
gives no contribution. On figure 1 we represent three possible situations. The fact 
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Figure 1: The quarter-plane of recurrent relations for D 

that the element of (2fc, q) is expressed in terms of other elements is reflected by a 
flesh. Let us note that the right-hand side of (3.3) uses also the moments M21 with 
21 <2k — 2. So, the flesh lines directed from points (2A;, ■) to {2k — 2, •) also indicate 
the dependence of D on M. 

Suppose we want to prove the estimate (3.6) for the elements D of the point 
{2k', q') with given k' and q' . On Figure 1 this is the point (10,5). We assume that 
these estimates are true for the elements D of all points {2k", q") that belong to the 
triangle 2k" + q" < 2k' + q' and of those situated on the line 2k" + q" = 2k' + q' with 
2k" < 2k' — 2. On Figure 1 we marked these points by fat circles. Also we assume 
that (3.5) holds for all M21 with 21 < 2k' - 2 + q' . Then we use relations (3.3) to 
deduce the estimate (3.6) for the point {2k', q'). 

So, starting with the point (2,2), we obtain the estimates (3.6) for the elements 
of (4,3) and (4,2). Then we go step by step over the points of B until we reach 
the point {2k', q') = (10,5). The next points are (12,4), (14,3) and so on. On this 
way we proceed along the line 2k" + q" = 2k' + q' until q" = 2 (the point (16,2) on 
Figure 1). Using (3.2), we prove the estimate (3.2) for M2k'+q'- Then we complete 
the next line starting from the point (10,6). 

Computationally this recurrent procedure means that we assume that all terms 
of the right-hand side of (3.3) verify inequalities (3.5) and (3.6) and then estimate 
their sum to show that it is less than the expression corresponding to the term of 
the left-hand side of (3.3). The same is true for (3.2) with the initial point of the 
recurrence Mq^'' = 1. We note also that M2^'' = v'^ that obviously verifies (3.5) 
with k = 1 because m2 = v"^. Let us have a look at D^^l that serves as the initial 
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point for (3.3). It is easy to see that 

^ e r ]_r^ jj Irj^ h\ = ^ (3.7) 

that certainly satisfies (3.6). 

Now let us consider (3.2). We want to show that if (3.5) and (3.6) are valid for 
the terms of the right-hand side of (3.2), then (3.5) is valid for M^^K The first term 
of the right-hand side of (3.2) is bounded by 



j=0 



/ 4(^-1)3^-^ ^ / 4k'\'-' 

+ — ^^^2 j ^ 2^^rn2k-2-2jm2j < \ l + j m2k- 

The second term of the right-hand side of (3.2) is bounded by 

\ / ai+a2=2fe— 2 \ / 

because 

2k-Z k-1 

J2{2k-2- j)j = J2{2i- If < {2k - 2f/2. 

3=1 i=l 

Gathering these two estimates, one obtains that 



Ak^\^~^ 



^2k — 



4(^-1)3 



m2k 



and (3.5) obviously follows. 

Now let us turn to (3.3). Let us denote 

We assume that (3.5) and (3.6) are valid for the all terms of the right-hand side of 
(3.4) and show that their sum is bounded by the expression standing at the right- 
hand side of (3.6). Let us consider two first terms of the right-hand side of (3.4). 
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Repeating calculations that lead to (3.7), we can write that the sum of the first two 
terms is estimated by 



ai— 3 oi— 3 

^ (ai - 2 - j) ruj m2k-2-j + ^ J "ioi-2-j m2k-ai+j 

3=1 j=l 



N1 \ N^J 

The terms in square brackets represent the first and the third parts of the sum 

[(ai - 3) ■ mim2k-3 + • • • + 2ma:,-4m2k-ai+2 + 1 • rna^-3m2k-ai+l\ + 

{mai-2m2k-ai + "T'ai-im2fc-ai-l + • • • + "l2fc-ai "T-ai -2} + 

[1 • m2k-ai+imai-3 + 2m2k-ai+2mai-4: + . . . + (ai - 3)m2jk-3mi] < 

2fc-3 

(ai - 3) X] "T,j'"T.2fc-2-j'- 

j'=l 

Certainly, all the sums run over even numbers, but we do not care about this because 
n^2j+i = and all relations are still true when regarding the sums over even and 
odd numbers. Thus we obtain that the sum of two first terms of the right-hand side 
of (3.4) is bounded by 

n(^;iV)(^l- A^mafc. (3.8) 

Now let us consider the third term of the right-hand side of (3.4). Assuming that 
(3.5) and (3.6) hold for corresponding M and D, we can write that 

^2''g%(<?+i)(j,ai-2-i,a2,...,a,)| < 

n(A,;iV)^m2fe_2 e'^^^^^^ < n(^,;iV)^m2fc. (3.9) 

The fourth term of the right-hand side of (3.4) is estimated by the same expression. 
The fifths term is estimated by 

Nl+l 2iV^ J V 27V^ J ^a^+a,-2m2k- 

1 ^ 1 

U{Ag;N)—v'^Y^ma,+ai-2m2k-ai-ai < U{Ag;N)—m2k. (3.10) 
ai ai 
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Finally, the last term of (3.4) is less or equal to 

2^ (ai + - 2) a2 • • • ag IH m2k-2 < 



U{A,; Nym2k-2 E ^^^^ < n(A; N)^^m2k. (3.11) 

1=2 



The sum of all expressions (3.6)-(3.9) gives us inequality 

\Dit.,aJ<m,;N) 



2 \ a? 2fe - 2 
1 +^ + 



Taking into account that ai < 2k — 2, it is easy to show that the sum of the terms 
in square brackets is strictly less than 1 provided A; < N^/^. This implies inequality 
(3.6). 

Lemma 3.1 is proved. 

3.3 Asymptotic behavior of the moments of GUE matrices 
3.3.1 Convergence of the maximal eigenvalue 

It is very well known the numbers m2k = ^12^(21;) determined by (3.12) given by 
Catalan numbers 

Elementary computations imply that 7^2^(21;) < (2'i;)^'^. Then it follows from Lemma 
3.1 that given e > 0, we have estimate 



(N) 
2k 



< (ivVlTey^ (3.13) 



for all k < N'^/"' provided Ak"'/N'^ < e. 

Taking into account inequality (/max)^'^ < Yld=i = Tri/^'^, we can write that 

Prob{W > 2v{l + £ } < ^\ < for ah k < N^'\ 3.14 

Regarding k S> log N and using the Borel-Cantelli lemma, it is easy to deduce from 
(3.14) that limsupjY^oo ^max < 2v with probability 1. This estimate together with 
the convergence (1.2) with probability 1 implies that Zmax — ^ 2v with probability 1 
as iV 00. 
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3.3.2 High moments and the scale at the spectral edge 

Let us compare our results with those already known for the moments of GUE with 
= 1/4. It is known that the moments M^"^ are given by the following recurrent 
relation P^ : 

^^{N) _ 2k-l (jv) 2fc - 1 2k -3 k{k-l) (;v) 

Explicit expression for ?7i2^ = m2A;(l) (3.12) implies that 

, _ 2A: - 1 , _ 2A; - 1 2A; - 3 , 

Then it is easy to deduce from (3.15) that M2^^ admit the following estimates (see, 
for example JH]): 

M(f < ( 1 + ^1 for all A., N. (3.16) 



Also one can write that 

1 + m',, for all A;, AT. (3.17) 

Inequality (3.16) means that the moments admit the power-like estimates by 

(2f)2'^(l + o(l)) in the limit 1 ^ k <^ jV^/^ and that this behavior can change on 
the regime 

Inequality (3.16) can be obtained by using the orthogonal polynomial approach 
(see the early paper by Bronk 6_, where the scaling at the edge of the semicircle 
distribution has been determined for ther first time); explicit asymptotic expressions 
were found in the seminal paper by Tracy and Widom It was shown that the 
fraction k^/N'^ is really the optimal one in the sense that one cannot decrease the 
exponent of k and increase that of N in (3.16). 

Basing on the relations (3.2) and (3.3), we did not manage to obtain estimates 
as precise as (3.16). Our result (3.5) implies that 

M2^f < 1^1 + m2k for k < N^/\ (3.18) 

Inequalities of this type are sufficiently powerful to estimate the maximal eigenvalue 
of H but they do not reflect the real scale of the eigenvalue distribution at the edge 
of the limiting spectrum. 
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Prom another hand, we did not used the orthogonal polynomial approach to 
obtain (3.2) and (3.4). Therefore the positive counterpart is that our approach is 
apphcable for more general ensembles of random matrices than the GUE. In the 
next section we show how our approach works in the case of band random matrices. 

4 Moments of band random matrices 

4.1 Main technical proposition and proof of Theorem 2.1 

Our main goal is to study the moments 




where are given by matrices (2.1). The first observation is that M^j^^^ = 

(see Lemma 5.2 of section 5). Thus we can consider the even moments M^''^^ only. 
The next observation is that in the case of band random matrices we cannot derive 
recurrent relations for M2k themselves as it was for the case of GUE. Instead we 
find a family of random variables that make a close system. Let us start with the 
random variables 

L,(.) = Lf.^) = [i7W)]L. 

Using integration by parts formula (see section 5), it is not hard to show that the 
mathematical expectation of L verifies the following identity 

EL2fc(x) = v^Y. ^{L2k-2-MLj[x]}, (4.2) 
i=o 

where we denoted 

s=l 

One can say that Lj[x\ represents a partial trace of normalized by h. So, one can 
expect that the variance of L goes to zero when 6 — > oo and that the mathematical 
expectation in the right-hand side of (4.2) factorizes. In what follows, we prove this 
factorization and estimate the variance of L[x\. 

Regarding mathematical expectation M2^'^^(x) = 'EiLx{k), we obtain equality 

Mir)(x) = .^E<'^l2,(-)<'^N+-^4l2(-), (4.3) 

J=0 
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where 

ai+a2=2k—2 

We recall that L° = L-BL. 

In the next subsection we will prove Lemma 4.1 that implies the following esti- 
mates that are true for all 2k < b^^^: 



sup M^^''\x)<(l + ^] m^k, (4.4) 



x=l,...,N 

where the family {r?7,2A;}A;eN determines the semicircle law (1.2). 
It follows from (4.4) that 

Using the estimate (3.14), we see that 

N 



Prob{ei)>2<l + £)}< 



(l + £)'= 



for all k < 6^/^ provided ik^/b^ < e. Choosing k such that k/logN — > oo, we obtain 
the bound 

limsup max ll'^^'^''] < 2v (4.6) 
with probability 1 provided b > (logA^)^. Theorem 2.1 follows form (4.6). 

4.2 Recurrent relations for the generalized moments 

Let us introduce random variables that serve as the elementary blocks for the closed 
system of recurrent relations. These are the products 

LR{a) (x) = . . . ^(y.-i)^a.^ ^ 

where R is given by two r-dimensional vectors A and Y; Ar = (ai, . . . ,ar), ai + 
. . . + ar = a, Yj. = (yi, . . . , y^), yi G {1, . . . , N}, and ^^^^ represents a diagonal 
AT-dimensional matrix [^^y^jst = 6st'ip{{s — y)/b). 
In section 5 (see Lemma 5.2) we prove that 

EL^(„)=0 if a = 2fc + l. (4.7) 
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In present section we prove the following proposition. 
Lemma 4.1 If 2k < b^/^, then 



< ELfl(„) <li + — \ rn2k if a = 2k, (4.8) 



and 



/n(^g;6)m2fc, if\Ag\ = 2k, 



if\Ag\ = 2k + l, 
(4.9) 

where 



and 

\A-q\ =ai + ... + aq. 



To prove estimates (4.8) and (4.9), we derive a system of recurrent relations 
that resemble very much relations (3.2), (3.3), and (4.3). Then we use the recurrent 
procedure described in subsection 3.2. So, we do not describe the details of the 
derivation but give the general description. 



Regarding I'fl(a), one can always assume that ai > 1. Therefore one can write 
that 




where we denoted *i = and 



/ _ f "1 - 2, if Z = 1 
"'"la,-l, if Z 7^1. 

Now, factorizing the mathematical expectation, one gets 

ELs(„)(x) =i;2^ E {(iJ°i-i*i...*,_i/f^)[x]}E {(iJ"'-i-%r ••*r-ii^"-) 

1=1 3=0 
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1=1 j=o 

ot'. — X \M 2 

If one accepts the symbolic denotation Y7i=\ S^j=o = Sl7=o ' then we can rewrite 
(4.9) in the form close to (4.2) 

a-2 a-2 

EL«(„)(x) = ^ E{Lj[x]}Y.{La-2-j{x)] + ^ E {L}[x]L°_2_j(x)}. 

J=0 J=0 

Now let us derive recurrent relations for D\ to simplify the formulas, we accept 
denotation D^'^\y,x) = E {Li{y)L2[x] ■ ■ ■ L°[x]} with obvious agreement that the 
vectors Ai,Yi that correspond to Lj are not necessarily the same for different i. 
Denoting T2,...^q = L2 • • • and mimicking computations of section 3, and those 
presented above, we obtain that 

ai— 3 

D^i){y,x) = v'' ^ ¥.{Lj[x\La,-2-jTl,J+ 
J=i 

(p) 1 

2 N <? P 1 / ^ 

^E^i'^^l^^EE E E{L^..L;_,(i^"^'''*...*i^"^^^-l-^■) X 

s,t=l p=2l=l j=0 ^ 

/ (1) 1 

iJJ*...^if"rp J L°+i---L°(i/"i ) . (4.10) 

J yt \ ^ sy 

After factorization of the mathematical expectation, the first term of the right-hand 
side of (4.9) produces four terms 

E {L}[x]L°[x] • • • L°[x]}E {La,-2-.j{v)} + E {L^,,_2-,/(y)^2N • • • i?N}E {Lj[x\}+ 

E{K^-2-j{y)L"j[x]L"2[x] ■ --mx]} - E{L^[x] • ■■mx]}E{Ll^_2-Aymx]}. 

(4.11) 

The estimates of corresponding sums v'^ YljLi are the same as those obtained in 
section 3 with the only difference that we take sup^ of the left-hand sides and 
replace in (3.8) and (3.9) the factor U{Aq;N) by Il{Ag;b). 

Let us consider the last term of (4.9). As one can see, it produces the terms of 
the forms ^MD^^"^) and ^D^*"^). Let us consider them in more details. For the 
first term we can write that 

See E E{L^..L^_lL^+,...L^}x 

p=2l=2 j=0 
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p=2 /=2 j=0 



62 



62 



mai+ap-2 < n(ylq;6)— m2fc. 



(4.12) 



Here we have used the fact that YJi=2 Sjio 1 = — 1- The remaining estimate is 



,,2 9 P "i 



p=2 Z=2 j=0 



p=2 ^ \ / 



-2 



(4.13) 



Now it is easy to gather the estimates of (4.11) with expressions (4.12) and (4.13) 
and obtain inequality 



sup|L>(«)(x,y)| <n(^,;6) 



1- — I +-7^ + 

ai 



a? 2A: - 2 



m2k- 



To complete the proof of Lemma 4.1, it remains to check that inequalities (4.8) and 
(4.9) are true for and D(2) (ggg (3 7^) 



5 Auxiliary propositions 

5.1 Integration by peirts formula for GUE 

Let us consider random matrices 

Hij = Oiij + i/3ij, 

where a and j3 are real independent gaussian random variables with zero mean value 
and variances 

Ea2. = t;V2, B0fj=v'/2. (5.1) 
We consider also the symmetric continuation A: 

~ \ aji, iii>j 
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and anti-symmetric continuation of B. Then H = A + iB is hermitian matrix with 
the probabihty distribution (1.6). Integration by parts formula imphes that 

E H^yG{H) =BAl^B ^ + lE B^^ E (5.2) 

Let us consider G{H) = {H^)st- Then we can write that 
dG I N 

j=i u,v=i 

Y: {{W-%x{H'-%t + {W)sy{H'-^)xt} . 

Similarly we get 

dG ' 



xy 

Substituting these relations into (5.2) and remembering (5.1), we obtain that 

E {H^y{H^)st} = E \H^y\^ E • (5-3) 

Regarding this formula, one can say that formally 

That is the way that we use the partial derivatives in sections 3 and 4. 



5.2 Even and odd moments 

To prove equality (4.7), we consider the mathematical expectation of the product 

E {h^s, . . . Hs^^t.^'^SHtot, ■ ■ ■ i^t.,«o4o^ • • • ^^0.1 • • • H,^r-} (5-4) 

for fixed values of Si,tj,... and odd sum of a^'s. It is clear that the number of 
factors H is odd and since wc have gaussian random variables, the mathematical 
expectation is equal to zero. The same reasoning shows that the left inequality of 
(4.8) is true. 
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